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CYLINDER MASS AIR FLOW PREDICTION MODEL 

FIELD OF THE INVENTION 
[0001] The present invention relates to mass air flow into an engine, 
and more particularly to an engine control system for estimating current 
mass air flow and for predicting future mass air flow into cylinders of an 
engine. 

BACKGROUND OF THE INVENTION 
[0002] The air to fuel (A/F) ratio in a combustion engine affects both 
engine emissions and performance. With current emissions standards 
for automobiles, it is necessary to accurately control the A/F ratio of the 
engine. Accurate control requires precise measurement and/or 
estimation of the mass air flow into the engine. 
[0003] Traditionally, engine air flow is measured with a mass air 
flow (MAF) sensor or calculated using a speed-density method. While 
MAF sensors are more accurate than speed-density calculation 
systems, they are also more expensive. An estimation-prediction 
method dynamically determines air flow into the engine using a 
mathematical model. While this method enables more precise A/F 
ratio control than traditional methods, inaccuracies may occur as a 
result of calibration difficulties. 

SUMMARY OF THE INVENTION 
[0004] Accordingly, the present invention provides a vehicle system 
to predict mass air flow into cylinders of an engine (CAF P ). The vehicle 
system includes a throttle position sensor that generates a current 
throttle position signal (TPS), a mass air flow (MAF) sensor that 
generates a current actual MAF into the engine signal, and a manifold 
air pressure (MAP) sensor that generates a current actual MAP signal. 
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A controller determines a current estimated mass air flow into cylinders 
signal (CAF E ), determines a MAF transient signal, and determines a 
MAP transient signal. The controller determines a CAF P signal based 
on the current CAF E signal, the current actual MAF signal, the current 
MAP signal, the current TPS signal, the MAF transient signal, and the 
MAP transient signal. 

[0005] In one feature, the MAF transient signal is based on a pre- 
defined MAF gain limit and the MAP transient signal is based on a pre- 
defined MAP gain limit. 

[0006] In another feature, the MAF transient signal is based on the 
current actual MAF signal and a prior actual MAF signal. The controller 
sets the MAF transient signal to zero if the MAF gain limit is greater 
than a difference between the current actual MAF signal and the prior 
actual MAF signal. If the MAF gain limit is less than a difference 
between the current actual MAF signal and the prior actual MAF signal, 
then the MAF transient signal is based on a difference between the 
current actual MAF signal, the prior actual MAF signal, and the MAF 
gain limit. 

[0007] In still another feature, the MAP transient signal is based on 
the current actual MAP signal and a prior actual MAP signal. The 
controller sets the MAP transient signal to zero if the MAP gain limit is 
greater than a difference between the current actual MAP signal and 
the prior actual MAP signal. If the MAP gain limit is less than a 
difference between the current actual MAP signal and the prior actual 
MAP signal, then the MAP transient signal is based on a difference 
between the current actual MAP signal, the prior actual MAP signal, 
and the MAP gain limit. 

[0008] In yet another feature, the controller schedules a select set of 
model coefficients based on a measured engine parameter. The 
controller determines the CAF P signal based on the select set of model 
coefficients. The select set of model coefficients is based on engine 
speed and MAP. 
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[0009] In still another feature, the controller determines the current 
CAF E signal based on a prior CAF P signal. 

[0010] Further areas of applicability of the current invention will 
become apparent from the detailed description provided hereinafter. It 
should be understood that the detailed description and specific 
examples, while indicating the preferred embodiment of the invention, 
are intended for purposes of illustration only and are not intended to 
limit the scope of the invention. 

BRIEF DESCRIPTION OF THE DRAWINGS 
[001 1] The current invention will become more fully understood from 
the detailed description and the accompanying drawings, wherein: 
[0012] Figure 1 is a functional block diagram of a vehicle including a 
controller that estimates current mass air flow and that predicts mass 
air flow (CAF P ) into engine cylinders; and 

[001 3] Figure 2 is a flowchart illustrating steps of a CAF estimation- 
prediction method according to the present invention. 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 
[0014] The following description of the preferred embodiment is 
merely exemplary in nature and is in no way intended to limit the 
invention, its application, or uses. For purposes of clarity, the same 
reference numbers will be used in the drawings to identify similar 
elements. 

[0015] Referring now to Figure 1 , a vehicle 10 is shown and 
includes an engine 12 and a controller 14. The engine 12 includes a 
cylinder 16 having a fuel injector 18 and a spark plug 20. Although a 
single cylinder 16 is shown, it will be appreciated that the engine 12 
typically includes multiple cylinders 16 with associated fuel injectors 18 
and spark plugs 20. For example, the engine 12 may include 4, 5, 6, 8, 
10, or 12 cylinders 16. 
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[0016] Air is drawn into an intake manifold 22 of the engine 12 
through an inlet 23. A throttle 24 regulates the air flow through the inlet 
23. Fuel and air are combined in the cylinder 16 and are ignited by the 
spark plug 20. The throttle 24 is actuated to control air flowing into the 
intake manifold 22. The controller 14 adjusts the flow of fuel through 
the fuel injector 18 based on the air flowing into the cylinder 16 to 
control the A/F ratio within the cylinder 16. 
[0017] The controller 14 communicates with an engine speed 
sensor 26, which generates an engine speed signal. The controller 14 
also communicates with mass air flow (MAF) and manifold absolute 
pressure (MAP) sensors 28 and 30, which generate MAF and MAP 
signals respectively. The controller 14 communicates with a throttle 
position sensor (TPS) 32, which generates a TPS signal. 
[0018] The controller 14 estimates current cylinder air flow (CAF E ) 
and predicts future cylinder air flow (CAF P ). Similar estimation- 
prediction systems are disclosed in commonly assigned U.S. Patent 
Nos. 5,270,935, issued December 14, 1993, and 5,394,331, issued 
February 28, 1995, which are incorporated herein by reference. The 
control system according to the present invention estimates cylinder air 
flow (CAF E ) into each cylinder. The controller 14 commands the fuel 
injector 18 for each cylinder based on CAF P to provide a desired A/F 
ratio within the cylinder 1 6. The controller 14 also may control ignition 
timing of the spark plug 20 based on the CAF E . 
[0019] The estimation-prediction system determines the CAF E 
based on prior predicted CAF's (CAF P ) and a current measured CAF 
(CAF M ). CAF m is preferably synthesized from other physical 
measurements such as MAP, MAF, TPS and RPM. It is anticipated, 
however, that a physical CAF sensor can be implemented to actually 
measure the current CAF. Calculation of CAF E is described in detail in 
U.S. Patent Nos. 5,270,935 and 5,349,331. 
[0020] Estimator correction coefficients are used in a weighted 
comparison. The estimator correction coefficients are pre-programmed 
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into memory and are predetermined in a test vehicle through a 
statistical optimization process such as Kalman filtering. The estimator 
correction coefficients are scheduled based on at least one engine 
parameter. Statistical optimization of the estimator correction 
coefficients provides that for a given engine operating point the 
estimator correction coefficients eventually achieve a steady state. As 
a result, the estimator correction coefficients may be determined off- 
line (e.g. in a test vehicle) and pre-programmed into memory. 
[0021] In accordance with the present invention, CAF P is 
determined based on the estimates, current engine parameters, a set 
of predictor coefficients, and transient behavior. Exemplary engine 
parameters include TPS, MAP, MAF, and engine speed (RPM). 
According to the present invention, the predicted CAF P is calculated as 
follows: 

CAF P (k+1) = a 1 CAF E (k) + a 2 MAF(k) + a 3 MAF(k-1) + b 1 MAP(k) + 
2 MAP(k-1) + b 3 MAP(k-2) + CiTPS(k) + c 2 TPS(k-1) + 
c 3 TPS(k-2) + d!UMAF(k) + d 2 UMAP(k) 

where k is the current time event, the component UMAF accounts for 
large MAF transients, and the component UMAP accounts for large 
MAP transients. To ensure steady-state accuracy, the predictor 
coefficients are constrained according to the following equations: 

ai + a 2 + a 3 = 1 
bi + b 2 + b 3 = 0 
Ci + c 2 + c 3 = 0 

The predictor coefficients di and d 2 are not constrained. The predictor 
coefficients are scheduled based on at least one engine parameter. 
For example, the controller 14 looks up the predictor coefficients within 
a particular schedule zone defined by RPM and MAP at time k. The 
predictor coefficients are difficult to calibrate in scheduled zones that 
feature a mix of small and large transients at steady-state. 
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[0022] To alleviate the difficulty of calibrating the predictor 
coefficients within the schedule zones, the components UMAF and 
UMAP are used. The component UMAF is governed by the following 
equations: 

UMAF(k) = MAF(k) - MAF(k-1) - MAFDEL 
if MAF(k) > MAF(k-1) + MAFDEL, otherwise 
UMAF(k) = 0 

where MAFDEL is a predetermined constant (gain limit) that 
differentiates between small and large transient behavior in MAF. If 
there is small transient behavior in MAF, then UMAF is set to zero. 
The component UMAP is governed by the following equations: 

UMAP(k) = MAP(k) - MAP(k-1) - MAPDEL 
if MAP(k) > MAP(k-1) + MAPDEL, otherwise 
UMAP(k) = 0 

where MAPDEL is a predetermined constant (gain limit) that 
differentiates between small and large transient behavior in MAP. If 
there is small transient behavior in MAP, then UMAP is set to zero. 
Thus, the components UMAF and UMAP enable accurate calibration of 
the predictor coefficients during small or large transient behavior. 
[0023] Referring now to Figure 2, the estimation-prediction control 
system will be described. The estimation-prediction control system 
determines a current CAF E based on a prior CAF P during an estimation 
loop. The engine 12 is operated based on CAF P and CAF E . A 
prediction loop determines CAF P for a future engine event based on 
the results of current engine operation. 

[0024] At step 100, control determines whether a CAF estimate 
interrupt is signaled. If false, control loops back. If true, control 
continues with step 102 and reads the current engine conditions (i.e. at 
time k) including TPS, MAP, MAF, and RPM. In step 104, the 
estimator correction coefficients are determined based on a MAP and 
RPM schedule, as described above. In step 106, CAF E (k) (i.e. current) 
is determined based on CAF P (k) and a weighted comparison of CAF 
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error (CAFERR). CAFERR is determined based on CAF P (k) and 
CAF M (k) and the estimator correction coefficients. 
[0025] In step 110, control enters the prediction loop by determining 
the predictor coefficients. The predictor coefficients are determined 
based on the schedule zones as described above. In step 112, control 
determines whether small or large transient behavior is occurring in 
MAF. If MAF(k) is less than or equal to the sum of MAF(k-1) and 
MAFDEL, small transient behavior is occurring and control continues 
with step 114. If MAF(k) is greater than the sum of MAF(k-1) and 
MAFDEL, large transient behavior is occurring and control continues 
with step 116. In step 114, UMAF(k) is set equal to zero. In step 1 16, 
UMAF(k) is set equal to the difference of MAF(k), MAF(k-1), and 
MAFDEL. 

[0026] Control continues with step 1 1 8 and determines whether 
small or large transient behavior is occurring in MAP. If MAP(k) is less 
than or equal to the sum of MAP(k-1) and MAPDEL, small transient 
behavior is occurring and control continues with step 120. If MAP(k) is 
greater than the sum of MAP(k-1) and MAPDEL, large transient 
behavior is occurring and control continues with step 122. In step 120, 
UMAP(k) is set equal to zero. In step 122, UMAP(k) is set equal to the 
difference of MAP(k), MAP(k-1), and MAPDEL. 
[0027] In steps 124 CAF P (k+1) is determined. CAF P (k+1) is used in 
a future estimation iteration to determine CAF E . Control exits the 
prediction loop and stores both calculated values and measured values 
in memory in step 128 for use in a future estimation-prediction iteration. 
In step 129, control operates the engine 12 based on CAF E (k) and 
CAF P (k +1) as determined in steps 106 and 124, respectively. In step 
130, the air estimate interrupt is cleared and control ends. 
[0028] Those skilled in the art can now appreciate from the 
foregoing description that the broad teachings of the current invention 
can be implemented in a variety of forms. Therefore, while this 
invention has been described in connection with particular examples 
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thereof, the true scope of the invention should not be so limited since 
other modifications will become apparent to the skilled practitioner 
upon a study of the drawings, the specification and the following 
claims. 



8 



